Long range Coulomb interaction in the ground state of bilayer graphene 
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We report on our studies of interacting electrons in bilayer graphene in a magnetic field. We 
demonstrate that the long range Coulomb interactions between electrons in this material are highly 
00 ' important. We show that in the unbiased bilayer (where both layers are at the same electrostatic 

■ potential) , the interactions can cause mixing of Landau levels in moderate magnetic fields. For the 

04 ■ biased bilayer (when the two layers are at different potentials) , we demonstrate that the interactions 



are responsible for a change in the total spin of the ground state for half-filled Landau levels in the 



. valence band. 



Monolayer graphene is a two-dimensional hexagonal crystal of carbon atonrs which exhibits a number of physical 
' and electronic properties that have made it the subject of intense study since it was first isolated in 2004 [1]. The 
existence of a gapless, conical low energy quasiparticle spectrum E = hvpk near the K points of the Brillouin zone 



^ ' (with the slope of the dispersion relation defined by the Fermi velocity vp w c/300) has lead to the possibility 

a ^ of realizing relativistic effects in a table-top sohd state experiment [2]. Additionally, the unusual chiral nature of 

"j^ the charge carriers was confirmed by the observation of an anomalous integer quantum Hall effect [sl with Landau 

a ' level (LL) spectrum _E„ — ^/nhvp / ^b, including a LL at zero energy and ^/B dependence of the spectrum on the 
I ■ 

'■^ I magnetic field (As is the magnetic length). From a more practical point of view, the unusually high electron mobility 

ch ■ n . 

Q ' and conductance of graphene at room temperature [J| makes this a very promising material for the fabrication of 

I electronic devices. 

I I Bilayer graphene (BLG) [5] consists of a pair of monolayers bonded by relatively weak dimer bonds perpendicular 

[ to the plane of the monolayer sheets (see Fig. [1]). In this material, both the conduction and valence bands have low 

lO ■ energy structure consisting of two quadratic branches seperated by the energy associated with the dimer bond, 71. 

r The lower conduction band and upper valence band are degenerate at the K points. The existence of chiral charge 
carriers with a Berry's phase of 2tt was confirmed in the observation of the integer quantum Hall effect Q where the 



low energy LL spectrum is app,„>,i„,.tely linear in the lield with ,. . ±r,..yH(S|T) tor „ > where is the 

cyclotron frequency. This spectrum includes a doubly degenerate LL at zero energy [7|. 

It has been predicted theoretically Q, 3, Q and observed experimentally (first by doping of the bilayer 1^1 , and 
then by applying electrostatic gates [1 1| ) that a gap can be induced in the low energy band structure of BLG by 
breaking the symmetry between the two layers. Switching off of the conduction current (and hence, proof of principal 
for switching devices made from so-called 'biased bilayer graphene') by sweeping the Fermi energy through the gapped 
region has been observed at low temperatures [12] , and this development has lead to a surge in interest in this system. 

The single particle theory of this material is well known 0, Qj S S| , it has been shown that the electron-electron 
interactions are significant in monolayer graphene [isl . [l^ . The Coulomb interaction has been studied in the unbiased 
(and hence, ungapped) bilayer {15^. while the biasing potential was considered in the context of a ferromagnetic 
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FIG. 1: The lattice structure of bilayer graphene. The upper (lower) lattice is shown by solid (dashed) lines. The Au (Bi) 
sublattices are shown with small, filled (large, open) dots; the B^-Ai dimers by large, filled dots, (a) The top-down view; (b) 
the side-on view projected between the two arrows in (a). 

transition due to short-range interactions in the mean-field approximation at zero magnetic field [l^. Also, certain 
collective modes leading to an intra-Landau level cyclotron resonance have been predicted in the presence of a magnetic 
field within Hartree-Fock theory 17]. However, the effect of the long range Coulomb interaction on the ground state 
of the biased system in a strong magnetic field has not been systematically investigated, and we address this problem 
in the current Letter. 

We find that the long range nature of the Coulomb interaction makes significant changes to the properties of the 
low energy charge carriers for BLG in a magnetic field. Specifically, the interactions are significantly stronger for 
electrons in the lowest LL, and this manifests itself in an observable way. It allows the possiblity of mixing of the 
LLs which were well-seperated in energy when the interactions were not considered. By calculating the explicit form 
of the ground state wave function, we show how this mixing fundamentally changes the nature of the ground state in 
the biased bilayer, and how this change manifests itself in the total spin of the ground state of the interacting system 
for half-filled LLs. 

In BLG, each sheet contains two inequivalent triangular sublattices of carbon atoms which we label Au, B^, Ai, 
and Bi for the A and B sublattices of the upper and lower layers respectively (see Fig. [T]). In the Bernal stacking 
arrangement, the inter-layer bonds consist of dimers formed from atomic orbitals associated with the Bu and Ai 
sublattices. The energy associated with this bond is denoted 71 and throughout this Letter we assume that 71 = 0.4eV. 
We allow for a static electric potential U to be applied between the upper and lower layers, so that the upper layer 
has potential J7/2, and the lower layer — J7/2. In a strong magnetic field we can write the tight-binding Hamiltonian 
using the four component basis {Au,Bi,Ai,Bu] with ^ = +1 in the K valley, and the basis {i?;, A^, v4;} with 
$ = -1 in the K' valley as 0] 
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where tt and tt^ are the operators corresponding to electron hops between neighbouring atoms (in opposite sublattices 
in the same layer) . The eigenvalues of Hamiltonian llj are found from the quartic polynomial equation 
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where n £ Z+. The energies t cx 71, (5 oc ^ and e| are measured in units of Additionally, we denote the band 
of a particular LL by placing a '+' ') after the level index for the conduction (valence) band. The wave functions 
associated with the LLs of this spectrum for n > I are given by 

where the functions (pn for 71 > are the magnetic oscillator functions in the Landau gauge and the coefficients are 
defined so that the overall wave function is normalized to unity. There are also levels with n = 0±, which have wave 
functions '00+ — e*'"'^ ((^0,0,0,0), and tpQ_ = e^^^ ((/Jq, 0,0,0) with e = ±5; and 0^0+ and are defined by Eqns. 
([2]) and ([3]) with the appropriate substitutions for n and ^. When [7 = 0, these four states are degenerate yielding 
the eight-fold degeneracy (including the factor of 2 for spin) seen in the integer quantum Hall effect in BLG Q]. 
Note that for a sufficiently strong inter-layer potential, the single particle energies described here can cross at certain 
magnetic field strengths. We include the fermionic properties of the electrons by constructing Slater determinants for 
the non- interacting many body state wave functions. 

To include the effects of the long-range Coulomb interaction we consider the Hamiltonian 

2 ^ e - r.j\ 

where the vectors rij label the positions of the electrons, and e is the dielectric constant of graphene. For graphene 
mounted on an Si02 substrate, esi02 = 2.5 [l^. Other methods of isolating graphene could be considered by changing 
the value of this constant. The effect of the interaction term is to mix the non-interacting states into some linear 
combination which minimizes the total energy of the system. Our analysis is conducted by employing the exact 
diagonalization scheme ^19j in which we calculate the linear combination of non-interacting basis states which gives 
the ground state of the Hamiltonian Ti = Ti.Q + Ticoui- This method entails dividing the infinite sheet into rectangular 
cells of dimension x Ly [20] and applying periodic boundary conditions to the wavefunctions at the edges of each 
cell. The matrix elements of the interaction over the single particle states described above are evaluated exactly by 
calculating the integrals over the spatial extent of the cell numerically. The interaction between the cells is taken into 



2l|. 



account by adding the Madelung energy of a charged lattice 

The number of single particle states included in the Hilbert space from which the non-interacting many body basis 
states are constructed is determined as follows. There are four relevant quantum numbers: The LL index n, the valley 
^, the spin and the momentum k = Tim/ L^- The values of the momentum are fixed when the boundary conditions 
are applied to the cell, and are labelled by < m < Af — 1 with M = L^Ly/ (27rA^). Which LLs we select is governed 
by the details of the system we wish to model, and M is set by computational restraints. 

In order to reduce the size of the many body system (and so to improve the speed of calculation), we observe that 
the Hamiltonian conserves the total momentum /i — X^iLi ''^^i niod M . Therefore, we can perform the diagonalization 
seperately for each value of /i which reduces the basis size to approximately the 1/Afth part. Finally, since there is no 
spin-dependent term in the Hamiltonian, the total projection of spin on the z axis Sz is also a good quantum number. 
We can fix Sz to its minimum value (which is for even N and i for odd N) whilst still being able to recover all 
eigenstates of the total spin operator S"^ [l4j| . 

We consider two different cases: Firstly we demonstrate the strength of the interaction by calculating the shift in 
the energy of each LL due to interactions for U — Q (an unbiased bilayer). Then we examine the system where the 
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-0.6443 
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Energy shift 
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-0.6316 


2 

-0.6222 


3 

-0.6148 


4 

-0.6085 



TABLE L Energy shift per electron due to the Coulomb interaction for integer fiUing factors in the 0± LL for (7 = 0. Energy 
units are e^/(eAs), the number of momentum states M — 3, and the magnetic field B — 3T. 



Landau level n 


1 + 


2+ 


3+ 


4+ 


Vn = \ 


-0.4766 


-0.5001 


-0.5242 


-0.5160 




-0.4705 


-0.4880 


-0.5176 


-0.5110 


I'n =3 


-0.4645 


-0.4759 


-0.5111 


-0.5061 


= 4 


-0.4584 


-0.4638 


-0.5046 


-0.5012 



TABLE II: Energy shift per electron due to the Coulomb interaction in LL n > 1 at fiUing factor Energy units are / (eAs), 
the number of momentum states M — 3, and the magnetic field B — 3T. 

filling factor is negative, the inter-layer potential is sizeable and the magnetic field strong. In this case, we observe 
changes in the total spin of the ground state as a function of U and the magnetic field strength B. 

We model the unbiased bilayer near half-filling by takng a single particle Hilbert space consisting of electrons in the 
0-1- and 0— LLs with all possible spin and valley states at ?7 = 0. Each integer value of the filling factor is simulated 
by taking the number of electrons N = {i' + 4)M, and we have M = 3. Table H] shows the results of diagonalizing 
the resulting many body Hamiltonian and evaluating the change in energy from the non-interacting ground state for 
integer filling factors. We notice that the energy shift per electron reduces slightly as the LL fills. 

In Table Ull we show the energy shift due to the Coulomb interaction for electrons in higher LLs (i.e. for levels 
with n > 1) in the conduction band. We have taken a single particle Hilbert space consisting of all spin and valley 
states within one LL. The filling factor Vn within LL n can range between (corresponding to an empty level) and 4 
(a filled level), so that Vn = 4 and i^n+i = describe the same overall filling factor. The number of electrons is set by 
N — VnM , and in order to allow direct comparison with the lowest LL we restrict ourselves to M ~ 3. Table Ull shows 
that the energy associated with the interaction of electrons is very similar in each of the higher LLs, and that the 
interaction energy per particle is slightly reduced as the LL is filled. We have verified that the results are identical in 
the valence band. 

Together, Tables HI and Ull show that the effect of the long range electron-electron interactions is considerable, and 
that for this value of the magnetic field {B = 3T) the shift in the higher LLs is only about two-thirds that of the 
lowest LL. In Fig. [2] we show the energy shift and absolute energy of filled LLs as a function of the magnetic field. 
The strength of the interaction scales with e'^/{e\B) oc ^/B with a roughly constant coefficient, while the LL spacing 
goes with ujc « B, so for lower values of the field, the n = 0± level crosses the 2— and 1— levels as shown in Fig. 
[UJb). The data shown here were calculated with e = 2.5, modelling graphene [igl deposited on an Si02 substrate. 
For suspended graphene (where e « 1), it is conceivable that the effect of the interaction would be even stronger. 
Additionally, the effect of the inter-layer potential is to bring together the valence band LLs with low index, so for 
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FIG. 2: (a) The energy shift per electron of filled LLs. (b) The absolute energy per electron of filled LLs showing the crossing 
between the n = 0± degenerate level and the higher LLs in the valence band. In both plots U — 0, and M — 5. 
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FIG. 3: The total spin of the ground state of the (a) u ^ -2 and (b) = -6 systems. M = 2 and N ^ {u + 8)M. The lines 
show the crossing points of the single particle states. The graining is due to the finite interval between data points. 



finite bias it is plausible that the interactions will cause significant mixing between these levels. 

Now we turn our attention to the system with an inter-layer potential applied, so that U ^ 0. With a finite gap 
between the 0-1- and 0— levels and non-zero filling factor, it is possible to consider the negatively-doped system by 
taking only those single particle states which are in the valence band. Therefore we select the single particle states from 
which to form the Slater determinants describing the non-interacting basis states by taking all spin and valley states 
of the 0— and 1— LLs. We have M — 2 and the number of electrons is related to the filling factor hy N = {ly + 8)M. 
Diagonalizing this system for half-filled LLs (so for v — —6 and v = —2), and calculating the expectation value of 
the total spin operator over the resulting ground state as a function of the magnetic field and the inter-layer bias 
gives the data shown in the plots in Fig. [3] |22||. We have also superimposed lines which represent the energy at 
which the single particle energy levels cross, as labelled in the legend. The plots show that there is an abrupt change 
in the total spin of the ground state, and a range of parameters where there is a non-zero polarization of the spin. 
This transition is not directly related to the crossings of the single particle states, since the position and slope of the 
transitions do not match the corresponding lines which we have superimposed on the plots. This effect is therefore 
purely due to the Coulomb interaction, and in particular to the exchange contribution which acts to minimize the 
energy of spin-polarized many body states. 
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FIG. 4: The occupancy of the single electron states in the intaracting many body ground state for each region of the plots in 
Fig. [3] The total z-projection of the spin is always zero since we fix this when forming the non-interacting many electron basis 
states. 

Figure |4] shows the occupation of the single particle levels in the interacting many body ground state. For simplicity, 
we display only the LL index of the states. The actual ground state is a coherent combination of several non- interacting 
basis states, where the combination of LL indices is consistent but different arrangements of momentum and valley 
states each come with identical prefactors in the linear combination. This figure shows that in the lower-right region of 
the parameter space considered, the electrons occupy as many of the 1— states as possible. As the parameters change 
to the upper-left region, the 0— levels become successively more populated. This demonstrates the fundamental effect 
of the inter-layer potential and the electron-electron interactions on the nature of the ground state. 

In conclusion, we have shown that the long range Coulomb interaction between electrons plays an important role 
near the Dirac point in BLG. In the unbiased case, the interactions will cause electrons in the doubly-degenerate 
lowest LL to be lower in energy than those in the 1— level for small-to- moderate magnetic fields. If an inter-layer bias 
is applied to split the valence and conduction bands, the electron-electron interactions precipitate a transition in the 
total spin of the ground state of half-filled LLs for certain ranges of parameters. This effect will have fundamental 
implications for the design of devices made from this material. 
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